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ABSTRACT 

We present new infrared observations of the emission/reflection nebula IC 405 obtained with the 
Spitzer Space Telescope. Infrared images in the four IRAC bands (3.6, 4.5, 5.8, and 8.0 [iva) and 
two MIPS bands (24 and 70 /im) are complemented by IRS spectroscopy (5 - 30 /im) of two nebular 
filaments. The IRAC (8.0 ^tm) and MIPS imaging shows evidence of a bow shock associated with the 
runaway 09.5V star, HD 34078, created by the interaction between the star and nebular material. 
The ratio of emission at 24 to 70 /im is higher in the immediate vicinity of HD 34078 than in the 
outer filaments, providing evidence for elevated dust temperatures (T^ > 90 K) in the shock region. 
The nebular imaging reveals that the morphology is band dependent, with varying contributions from 
aromatic emission features, H2, and dust emission. Nebular spectroscopy is used to quantify these 
contributions, showing several aromatic emission bands between 6-14 /im, the S(5), S(3), S(2), and 
S(l) pure rotational emission lines of H2, and atomic fine structure lines of Ne, S, and Ar. The 
low-dispersion spectra provide constraints on the ionization state of the large molecules responsible 
for the aromatic infrared features. H2 rotational temperatures of the two bright nebular filaments are 
determined from the observed line strengths. An average T(H2) ~ 400 K is inferred, with evidence 
for additional non-uniform excitation by UV photons in the intense radiation field of HD 34078. 
The photoexcitation hypothesis is supported by direct measurement of the far-UV H2 fluorescence 
spectrum, obtained with FUSE. 

Subject headings: dust — ISM: molecules — ISM: individual (IC 405) — reflection nebulae — Infrared: 
ISM 



1. INTRODUCTION 

IC 405 is an emis sion /reflection nebula in Au- 
riga l|Herbidll958l 11999ft . It is illuminated by the 09.5V 
star HD 34078 (AE Aur, V = 6.0), a high proper mo- 
tion (« 40 mas/yr; Blaauw & Morgan 1954) runaway 
star ejected from the Orion star-forming region in a 
binar y-binary interaction ap proximately 2.5 million years 
ago ijBagnuolo et alJfeOOH) . The star an d the optical 
nebula are roughly cospatial at pres ent l)Herbid [l958t 
iFrance et all 120041 fBoisse et all 1200.51) . at a distance of 
450 pc. The sightline to HD 34078 is considerably red- 
dened (E{B -V) = 0.53; Cardelli et al. 1989) , however, 
the extended nebul ar regions appear t o be free of fore- 
ground attenuation l|France et alJl2004ft . 

The large transverse velocity of the star indicates that 
it has only recently encountered the nebular material 
with which it now interacts. The present geometric 
configuration has made the HD 34078 sightline a tar- 
get for studie s of interstellar absorption line variation. 
iHerbid l)1999ft combines fifty years of atomic and molec- 
ular line observations in the violet-blue (A3400 - 4900 
A) to set an upper limit on the length scale of sub- 
structure in the ISM of 1450 AU. More recently, vari- 
ations in the column density of interstellar CH and CN 

1 Based in part on observations made with the the 
Spitzer Space Telescope, which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under a contract 
with NASA. 

2 Current Address: Canadian Institute for Theoretical Astro- 
physics, University of Toronto, 60 St. George Street, Toronto, ON 
M5S 3H8 



(10 - 20 %) have been observed on this sightline, set- 
ting a smaller upper limit on the substructure siz e (< 
200 AU) tow ards HD 34078 llRollinde et all 12003ft . A 
similar study ijBoisse et al.l f2005) using Far Ultraviolet 
Spectroscopic Explorer (FUSE) observations has set 
an upper limit on the molecular hydrogen (H2) column 
variation of < 5%, constraining the scale to a few tens of 
AU. They also find evidence for interaction between HD 
34078 and the surrounding nebula. Detections of excited 
molecular hydrogen absorption lines (measurements up 
tov = 0,J=ll and tentative detections up to the v = 2 
level) are evide nce for fluorescent pumping by the local 
stellar UV field {B oisse et aT1l2005ft . Additional observa- 
tions of the interaction between the UV radiation field 
of HD 34078 and the ne bular gas and dust are presented 
by IFrance et alJ {2004) . They present a spectroscopic 
study of far-UV dust scattering and fluorescent H2 emis- 
sion lines in the nebula. 
Evidence for interaction between the star and 

nebular dust grains h as also been detected by 

IRAS Ijvan Buren et all 11995ft . The high proper 
motion of HD 34078 makes it a candidate host for a 
bow shock. As a massive star moves supersonically 
through an ambient interstellar medium, stellar radia- 
tion pressure and/or strong stellar winds can produce 
a comet-shaped, compact H II region surrounded 
by s hocked molecula r gas ijvan Buren fc McCravl 
119881: Ivan Buren et all 11990ft A 60 itm excess 
observed in IRAS map s l)van Buren et alJ Il995t 
iNoriega-Crespo et all 11997ft identifies IC 405 as a 
potential bow shock region. 
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In this paper, we present the first direct imaging of the 
bow shock around HD 34078 with the Spitzer Space 
Telescope-IRAC and MIPS. We also use MIPS data 
to study the nebular dust emission, characterize the 
spectroscopic properties of aromatic emission features 
(AEFs) and H2 using nebular observations made with 
the Spitzer-IRS, and put these new observations in con- 
text with existing studies of IC 405. §1.1 provides back- 
ground information on the diagnostic properties of AEFs, 
H2, and dust in diffuse nebulae. In §2, we describe the 
new Spitzer data, and introduce supporting far-UV spec- 
troscopy from FU SE. An analysis of the infrared (IR) 
imaging and spectroscopy is given in §3. In §4, we dis- 
cuss the structure of the bow shock region, present de- 
terminations of the H2 rotational and dust temperatures, 
and analyze spatial variations in AEF carrier ionization 
across the nebular filaments. Our findings and a sum- 
mary of the observations are presented in §5. 

1.1. AEFs, Hi, and Dust Emission in Photodissociation 

Regions 

AEFs have been found to be ubiquitous in the mid-IR 
spectra (3 - 20 fim) of reflection nebulae and H II regions 
(see van Dishoeck 2004 for a recent review). AEFs are 
seen as broad bands in the /SO-SWS/JS'O-CAM-CVF 
and Spitzer-IRS wavelength ranges, with lines com- 
plexes centered on 3.3, 6.2, 7.7, 8.6, 11.2, and 12.7 /im. 
These features are usually attributed to a superposition 
of C - C and C - H vibrational modes in a distribu- 
tion of polycyclic aromatic hydrocar bons (PAHs) consist- 
ing of tens to hund reds of C-atoms l)Schutte et alJll993t 
Ivan Disho eck 2004). However, discrepancies between as- 
tronomical emission features and laboratory PAH spectra 
exist, preventing a conclusive identification of the exact 
carrier of the AEFs (Gordon et al. 2006 and references 
therein). The carrier population that gives rise to the 
observed mid-IR bands (also called Aromatic Infrared 
Bands; Verstraete et al., 2001) spans a range of sizes 
and ionization states. The size distribution is thought 
to be regulat ed by the hardness of the illuminating ra- 
diation field ijVerstraete et aljl2001f) while the aromatic 
molecule ionization balance is controlled by the strength 
of the u ltraviolet (UV) radiation field and the el ectron 
density l|Bakes et al.ll200l iBreginan fc Temil ETOQ.^. The 
relative strengths and spectral profiles of the major AEFs 
reflect the physical state of the carrier population in a 
given environment. 

Variations in the spectral features of the AEFs 
can be used to quantify the response of the carrier 
molecu les to different radiation and density environ- 
ments (|Bakes et al.ll200l|) . These variat ions can be stud- 
ied b y comparing a sample of nebulae l) Verstraete et al.l 
120011). at different locations within a single ob- 
iect l)Verstraete et al.lll996t [ Honv et al.ll2001|). or in the 
diffuse interstellar medium l lB^ukmlreret"liL| 2000'). The 
C - H emission features at 3.3 and 11.2 /jm are found 
to be well correlated over a range of physical conditions, 
however the relative strengths of the C - H an d C - C fea- 
tures are observed to vary with environment l)Honv et al.1 
l200ltlvan Dishoeckl2004|) . The ratio of a C - C to a C - H 
feature (6.2/11.2 /im for example) reflects the ionization 
state as C - C band strengths are stronger in ionized 
aromatic carri ers while C - H ba nds display the oppo- 
site behavior ljUchida et al.ll2000j) . The increase in the 



6.2 /zm C - C feature (and more generally, the 6.2, 7.7, 
and 8.6 /iia bands) is due to the dipole moment pro- 
duce d by a change in ch arge distribution during ioniza- 
tion ijPeeters et al.ll2002|) . The carriers of the AEFs fall 
in the middle of the molecular mass distribution in reflec- 
tion nebulae and H II regions (which we will group under 
the broad heading of photodissociation regions; PDRs). 

H2, residing at the lower limit of this distribution, 
make s up the majority of the molecular mass in thes e ob- 
jects ijHollenbach fe TielensllT997l iShaw et al.ll2005ft . H 2 
is seen to emit from far-UV to the mid-IR wavelengths, 
excited by both UV photons (i.e. - fluorescence) and 
thermal processes, such as collisions and shocks. Detec- 
tion of the UV emission lines of H2 is one clear indication 
of photoexcitation. H 2 dissociates in environments with 
temperatures corresponding to the exc ited electronic lev- 
els th at give rise to the UV lines ijShull &: Beckwithl 
1982). The emission nebula IC 63 is the prototypical ob- 
ject for PDR studies of fluore scent H 2 . The far-UV spec- 
trum has been well-stud ied l)Witt et al.l Il98flt iHurwitzl 
119981 iFrance et al.N2005|) . confirming stellar continuum 
photons as the excitation source of the observed emis- 
sion. F ar-UV H? emission lines in IC 405 have been ob- 
served l|France et al. 1 12004]) . and these observations will 
be discussed below. 

Near-IR observations of the rovibrational emission 
lines of H2 can also be used to constrain the molec- 
ular excitation mechanism. Collisional processes can 
populate the lowest excited vibrational states of the 
molecule (v — 1 and 2), but the fluorescent cascade 
associated with the electronic excitation of H2 by UV 
photons is the only mechanism (ignoring excitation by 
non-thermal electrons) available to populate higher vi- 
brational states (v > 2) as the molecules dissociate 
at kinetic temper atures proportional to these higher 
vibrational levels l|Shull fc B cckwith 11982ft . Ratios of 
these excited vibrational states can be compared with 
model predictions to determine th e excitation mecha- 
nism (Bl ack fc van Dishoeck! Il987ft . N ear-IR line ra- 
tios indicate UV excitation in IC 63 fL uhman et alJ 
IT997]) the bright reflection nebulae NGC 2023 and NGC 
7023 l|Martini et alJ 119991: iTakami et all 120001) . and o n 
larger scales in molecular clouds (|Luhman fc .T affc 1996). 
Collisions contribute to the rov ibrational excitation o f 
H 2 in dense molecular clouds llKristensen et al.1 12003|) 
and regions where sho cks play an importa n t role, such 
as pl anetary nebulae (K astner et aljTl996T: ILupu et all 
2006). The pure rotational emission lines of H2, ac- 
cessible only in the mid-IR, give diagnostic information 
about the molecular gas temperature in PDRs. Ra- 
tios of these lines have been u sed to derive temper- 
atures in star-forming re gions (lRo s^nthalet_alJ |2000|; 
Alters et al.ll2005|) . PDRs llThi et al.lll999UHabart et alJ 



and external gal axies IjValentiin fc van der W cri 



1999HArmus et aDl200l . 

Dust grains populate the upper end of the molecu- 
lar mass distribution in PDRs. Grains are observed 
through their thermal emission in the mid/far-IR (Schnee 
et al. 2005 and references therein) and in the atten- 
uatio n and scattering of ultraviolet and visible pho- 
tons dCardelli et al.lll989t iWitt et aTlll993t iBuYgh et alJ 
120021 iDrainei 12003ft . IRAS observations of PDRs find 
that grain e mission is spatially cor r elated with near-IR 
H 2 emission l)Luhman fc ,Iaffell996T) . I.Iansen et all (|1994l 
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Fig. 1.— Left: Optical image of IC 405 from lFrance et al.l<200l) . Right: Post-BCD 8.0 /im IRAC images of IC 405, scaled to highlight 
the rich nebular structure. The origin is the position of the illuminating, runaway 09.5 star, HD 34078. The main nebular regions studied 
in this work are labeled A and B. White overlays represent far-UV spectrograph apertures (the squares are the FUSE LWRS slit and the 
rectangles represent a rocket-borne spectrograph; France et al. 2004). At the cores of the bright filaments, discrete features contribute < 
40% of the nebular emission at 8 /im. 



1995) present observations of IC 63, where dust emission 
is coincident with the region of peak H2 emission (at UV 
and IR wavelengths) as well as CO and other molecules. 
This correlation would be expected as dust grains are 
thought to be the p rimary site of H2 formation in the 
interstellar medium l|Cazaux fc Tielensll2002t 120041) . Ad- 
ditionally, photoelectric emission from small grains is a n 
important source of heating in PDRs ijAbel et al.ll2005j) . 
lHabart et alJ l)2004D provide observational evidence for a 
correlation between AEFs and H2, which they propose 
is related to H2 formation on PAHs. With the wealth of 
observational evidence that these molecules are interre- 
lated in PDRs, it follows naturally to study this group 
as a whole. 

In this work, we further address the physical conditions 
of the AEF carriers and their relationship with H2 and 
dust grains by presenting new observations of the emis- 
sion/reflection nebula IC 405 made with Spitzer. We use 
a combination of imaging and spectroscopy to constrain 
the variation in aromatic carrier properties with chang- 
ing radiation and dust environment through the nebula. 
Our aim is to present new IR observational characteris- 
tics of large molecules, dust grains, and H2 in the unique 
PDR environment of IC 405, for comparison with previ- 
ous far-UV observations. 

2. OBSERVATIONS 
2.1. Spitzer 

Spitzer observations of IC 405 (Program ID 20434) 
were made at positions coincident with the UV data pre- 



sented in lFrance et al.l {2004), where H2 fluorescence was 
observed above the level of dust scattered light (FU SE 
Pos3: RA = 5 h 16 m 30.21 s , S = +34°24'56.4" and FUSE 
Pos4: RA = 5 h 16 m 18.64 s , S = +34°32'25.2", J2000; see 
also §2.2). The Infrared Array Camera (IRAC) observed 
the nebula in mapping mode on 2005 September 23, with 
an exposure time of 12 seconds per frame, for a total 
integration time approximately 825 seconds. Details o f 
the IRAC instrument can be found in lFazio et alJ {2004) . 
The IRAC pointings covered the northern nebula in all 
four bands (rl6033792), with the 4.5/8.0 am channel ac- 
quiring data of HD 34078 and the associated bow shock 
(rl6033280, Figure 1). The images were reduced with 
the data processing pipeline (S14.0.0) provided by the 
Spitzer Science Center (SSC). Analysis was performed 
on the post-Basic Calibrated Data (BCD) products. 

Using the mapping mode on the Infrared Spectrograph 
(IRS, Houck et al. 2004), we obtained spectra clustered 
around previous FU SE observations Pos3 (rl6033024) 
and Pos4 (rl6033536). IRS programs were carried out 
on 2006 March 17 and 19 with a total exposure time 
of ~ 9600 seconds. Data at all IRS wavelengths were 
acquired through both the high and low-resolution aper- 
tures (SL, SH, LL, and LH). The spatial coverage of the 
low-resolution aperture scan allows us to correlate spec- 
tral variations with nebular morphology, hence this work 
will focus on an analysis of the low-resolution (Short- 
Low, 5 - 15/im; R = 60 - 130) data. Short-High and 
Long-Low spectroscopic data were used as supporting 
observations. The post-BCD (pipeline version S13.2.0) 
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TABLE 1 

Mid and far-IR. photometry of bow shock region. 



Position 


RA (2000) a 


8 (2000) 


F 8 


F 24 


F 7 o 






^ h m sj 


(° ' ") 


(Jy) 


(Jy) 


(Jy) 


(K) 


$E 


05 16 20.98 


+34 18 52 


0.508 ± 0.083 


3.435 ± 0.229 


1.945 ± 0.169 


83.93 ± 2.06 


S C 


05 16 18.42 


+34 18 58 


0.805 ± 0.103 


4.286 ± 0.255 


1.904 ± 0.167 


88.76 ± 2.24 


s w 


05 16 16.83 


+34 18 48 


0.675 ± 0.094 


3.347 ± 0.225 


1.805 ± 0.162 


84.85 ± 2.16 



a Shock pointings were defined from analysis of the 24 /im MIPS data obtained near HD 34078. 
b Dust temperatures were derived from the ratio of 24 to 70 fim flux density, F24/F70. This 
temperature should be considered a lower limit due to the saturation of the 24 /jm map. 



reduced one-dimensional spectra were tagged according 
to the central aperture position for registration with the 
associated filaments observed in the imaging. As the "on- 
filament" nebular positions filled most of the Short-Low 
aperture, analyzing the full-slit extraction at all point- 
ings yielded information about both the spectral charac- 
teristics and the slit-filling fraction. 

We also obtained mosaic maps of IC 405 in the three 
bands of the Multiband Imaging Photometer for Spitzer 
( MIPS). A des c riptio n of the MIPS instrument is given 
in lRieke et all pOO^ . The data (rl6032768) were ob- 
tained on 2006 February 28. The post-BCD scan maps 
(S13.2.0) at 24 and 70 /im (Figures 2 and 3) were used 
to perform photometry of several nebular filaments and 
complement the imaging at shorter wavelengths. The 24 
and 70 pirn maps are of high quality, but we note the 
inner 10 - 25" of the bow shock region near HD 34078 is 
likely saturated on the 24 /im Si:As array. The 160 /im 
map contains gaps that make measuring long- wavelength 
photometry very challenging, we will focus on the shorter 
wavelength MIPS channels here. The 24 and 70 /im scans 
covers roughly 8' x 54' and 5' x 52', respectively, both 
centered on FUSE Pos4. 

2.2. FUSE 

IC 405 was observed by FUSE on 2003 March 11-13 
as pa rt of the D127 guest observer program ijFrance et alJ 
,2004). Spectroscopy was obtained in the FUSE far- 
UV bandpass (905 - 1187 A) through the low-resolution 
(30" x 30"; AA w 0.33A) aperture (Moos et al. 
2000). The data acquired at positions Pos3 and Pos4 
show H2 emission lines and were used to guide the 
Spitzer observations presented in this work. We use the 
LiF2a (1087 - 1179 A) channel to measure the strength 
of the strongest far-UV band of continuum pumped H2 
fluorescence, A ~ llOOA, for comparison with the pure 
rotational emission lines in the IRS bandpass (§4.2). The 
FUSE data were obtained in time-tagged mode and re- 
duced with the CalFUSE pipeline, version 2.2.3, with 
errors determined from the associated count rate plots. 
Both positions were observed for roughly 12.5 ksec. 

3. ANALYSIS AND RESULTS: IMAGING AND 
SPECTROSCOPY 

The 8.0 /im IRAC images showed the best combination 
of S/N (~6 - 38 for the faintest to the brightest nebular 
features), spatial resolution (<2"; Fazio et al. 2004), and 
nebular structure with a field-of-view that covered all of 



TABLE 2 

8 (IM IRAC PHOTOMETRY OF NEBULAR FILAMENTS. 



Position 


Program 


RA (2000) 


5 (2000) 


Flux a 






^ h m 


(° ' ") 


(Jy) 


A 


rl6033280 


05 16 30.97 


+34 25 04 


0.466 ± 0.079 


A* 


rl6033280 


05 16 19.84 


+34 25 18 


0.385 ± 0.072 


A** 


rl6033280 


05 16 15.53 


+34 23 57 


0.573 ± 0.088 


B 


rl6033792 


05 16 16.88 


+34 32 21 


0.543 ± 0.085 



a The 8.0 /mi images contain contributions from AEFs, H2, 
and [Ar II]. The 7.7 AEF contributes between 30 - 40% of the 
8.0 /im band flux in the region of peak filament brightness. 



the subregions of interest. A composite of the 8.0 fim 
images is shown in Figure 1. These data provided a refer- 
ence frame for all of the image and spectroscopic analysis 
described below. The overall morphology of the 8.0 /im 
images shows the brightest structures seen in the three- 
color optical image 3 , displayed in Figure 1 for comparison 
with the IRAC data. The "shell" structure seen to the 
East/Northeast of HD 34078 with a radiu s of ~ 400" ob- 
served in DSS images and in Plate IV of iHerbi el l|1958f) 
is observed at 8.0 /im. 

Two bright nebular regions were chosen for analysis in 
this work, based on previous ultraviolet studies. Fila- 
ments A and B were part of the [France et alJ l)2004|) far- 
UV study of IC 405. Data acquired at these positions 
showed a combination of fluorescent H2 emission and 
dust-scattered starlight. These observations provided de- 
tection of continuum pumped H2 emission where the ro- 
tational substructure was spectroscopically resolved at 
wavelengths shortward of Lya. Filament A is broken 
up into three subregions: the main A knot with exist- 
ing UV spectroscopy; ~175" to the west, Filament A* 
is a clumpy feature located at the end of the Filament 
A structure; and Filament A** to the southwest was the 
brightest feature in the 70 fxm data (these nebular fea- 
tures are labeled in Figure 1). The brightest features in 
the optical image are related to the regions containing 
Filaments A and B at infrared wavelengths. Filaments 
A and B are also the dominant features in the 6300 - 

3 T.A. Rector, http://antwrp.gsfc.nasa.gov/apod/ap011204.html 
also presented in l igure 1 ot trance et al. (^0U4). 
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Fig. 2. — Post-BCD 24 /jm map of IC 405. The region near 
HD 34078 (0,0) shows the wake-like morphology, suggesting the 
presence of a bow shock as the star traverses the ambient nebu- 
lar material. The bow shock region at 8 and 24 fim can be seen 
enlarged in Figures 4 and 5. The overlays on Filament A and B 
represent the IRS-SL scan regions, with the ON/OFF convention 
used in Figures 10 and 12 indicated. (The IRS-SL scans are shown 
in more detail in Figures 8 and 9.) 



Fig. 3.— Post-BCD 70 /an scan of IC 405. The position of HD 
34078 is at the origin. The Filament A** region is the strongest 
feature in the 70 fim data. 



6750 A photographic images presented in Herbig (1958; 
Plate III), where line emission from Ha, [N II], and [S II] 
is the primary source of the nebular brightness. 

The following subsections are divided into imaging and 
spectroscopic results. The imaging is presented first, 
starting near HD 34078 and moving out into the north- 
ern filaments. An analysis of the IRS spectroscopy fol- 



6 



France et al. 



4-00 



300 ■ 



200 



0: 



-100 





A" 






IRAC- 8.0 |im 








"Cirrus" 













-400 



-300 -200 -100 
Arc seconds 



Fig. 4. — IRAC 8.0 /im blow up of the bow shock region, scaled 
to illustrate the shock structure. The shock subregions are labeled 
Shock-East (S B ), Shock-Center (So), and Shock-West (Sw). Pho- 
tometry of the subregions is given in Table 1. 
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Fig. 6. — IRAC 4.5 /im image of the HD 34078 region with a 
logarithmic intensity scale to enhance the nebular emission. The 
"cirrus" features described in §3.1 are identified. The 4.5 /an band 
shows these structures at a higher contrast than the shell structure 
seen in Figure 1. 
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Fig. 5. — MIPS 24 (im blow up of the shock region, scaled 
to illustrate the shock structure. The stellar position is at the 
origin. The inner 10 - 25" near HD 34078 in the direction of stellar 
motion (North, +y) is saturated at 24 fim, indicative of the high 
temperature of the grains in the shock region. The dashed lines 
indicate the location and direction of the spatial profiles presented 
in Figure 7. 24 fim photometry of the shock region is given in 
Table 1 



lows, we determine the relative contribution of the neb- 
ular species to the IRAC images and observe position- 
dependent line strengths. We note that the IRAC band 
overlap regions and spectroscopic programs were cen- 
tered on Filaments A and B, thus six-band imaging and 
spectroscopy do not exist for all IC 405 subregions. A 
brief comparison between the imaging and spectroscopy 
of the different regions is given in §3.3. 

3.1. Imaging: IRAC and MIPS 

Bow Shock - The wake feature seen clearly in the 8.0 
/im IRAC (Figure 4) and 24 /im MIPS (Figure 5) im- 
ages is thought to be enhanced dust emission produced 
by the in teraction of HD 34078 with t he ambient nebular 
medium Jvan Burc nfc McCravM l988'l. The bullet-tip in- 
teraction region is seen most clearly in the 24 /im images 
(despite the saturation of the brightest ridge) , and three 
positions were chosen for analysis. Moving from left to 
right on the image blow-ups, they are labeled Shock-East 
(S B ), Shock-Center (S c ), and Shock- West (S w ). The 8.0 
/zm data show an emission arc of diameter ~ 75" with 
an enhancement of factors of > 5 (in units of MJy sr _1 ) 
relative to the interfilament medium. The emission arc 
has a brightness peak 10 - 20" to the northeast of the 
stellar coordinates. The 24 /im image shows a similar 
75" arc, with higher contrast emission in the outer edges 
of the arc. The 24 /im data also suggest a connection 
between the shocked material and the cirrus features dis- 
cussed below. The 24 /im flux density transitions from 
the cometary arc into the cirrus features to the east of 
the star more smoothly than at 8 /im. Both the 8 and 24 
fim images show a clear spatial offset between the posi- 
tion of HD 34078 and the peak of the infrared emission. 
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Fig. 7. — Spatial profiles of the 24 and 70 /im MIPS scans of the 
southern region of IC 405. The 24 fim data have been convolved 
to a resolution of Rj 18" for comparison with the 70 fim map. The 
direction and magnitude of the stellar motion through IC 405 is 
indicated. The 24 (im profiles outline the shock region while the 
70 /im cuts show a weaker peak near the stellar position and a 
maximum near the Filament A** region. 



TABLE 3 

24 fJ,M MIPS PHOTOMETRY OF NEBULAR FILAMENTS. 



Position Program RA (2000) <5 (2000) Flux a 

( h m s ) (° ' ") (Jy) 



A rl6032768 05 16 30.99 +34 25 03 1.045 ± 0.126 

A* rl6032768 05 16 19.70 +34 25 18 0.984 ± 0.123 

A" rl6032768 05 16 15.52 +34 23 58 1.175 ± 0.134 

B rl6032768 05 16 18.17 +34 32 27 0.892 ± 0.117 



a IRS Long-Low spectra were analyzed to assess the emission 
line contribution to the 24 fim band flux. Discrete lines were 
found to contribute less than 8% of emission in the 21.5 - 26.2 
fim bandpass. 



The 4.5 /im IRAC image is dominated by the saturated 
core of HD 34078 (Figure 6), but asymmetries in the 
PSF suggest that an enhancement exists in the direction 
of stellar motion. An arc is detected in the 70 /im MIPS 
observations, which can be seen faintly in Figure 3. 

We measured the variation in dust profile across the 
interaction region by taking spatial cuts through the 24 
and 70 /im MIPS data, intersecting the Sg, Sc, and Sw 
positions. In order to make a fair comparison, the 24 /im 
data were convolved to a resolution of 18" and rebinned 
to the dimensions of the 70 /im scan. The images were 
then aligned with a ~ 3° rotation to account for space- 
craft roll orientation and registered with their WCS data 
using a combination of image analysis tools (custom IDL 
codes as well as ds9 from the Chandra X-Ray Center and 
the SSC data retrieval program Leopard). We found that 
the registration was accurate to roughly 2" in x (~ RA) 
and 10" in y (~ Dec). These profiles are shown in Figure 
7 running from south to north (left to right on the figure), 
intersecting the Se, Sc, and S^ coordinates listed in Ta- 



ble 1. Se, Sc, and Sw are labeled in Figures 4 and 5, 
with dashed lines representing the cuts in Figure 5. The 
24 /im profiles show the dominance of the bow shock, 
again noting that the peak brightness regions (10 - 25" 
from HD 34078) are likely saturated. The position of HD 
34078 is indicated with the dotted line at and the direc- 
tion of stellar motion through the nebula is labeled in the 
center panel. The Se profile shows the largest breadth 
due to its association with the cirrus-like structures. The 
70 /im profiles show an enhancement at the shocked re- 
gion, but rise to a maximum at the intersection of the 
nebular arm containing Filament A** (y « 400", marked 
with the dashed line in Figure 7). The second peak in 
the 70 /im Se and Sc profiles near 520" are associated 
with the intersection of Filament A*. We also observe 
an additional peak at the intersection of the Filament 
B arm, but we do not show a comparison as the image 
registration was of poorer quality. 

Flux ratio (F24/F70) profiles were made using these ex- 
tractions. These profiles were interpolated onto a com- 
mon spatial grid (« 4.7"/pixcl) before comparison. The 
registration in y was the primary source of uncertainty 
for these comparisons. These ratios were used to derive 
the dust temperature profile in the shocked region, dis- 
cussed in §4.1.1. Photometry of the Se, Sc, and Sw 
regions was performed on the 8, 24, and 70 /an images 
following the procedure described below. The photomet- 
ric values for the shock regions are listed in Table 1. 

The HD 34078 Region - The largest difference be- 
tween the 8.0 am and optical images of IC 405 are the 
cirrus-like sheets of emission that are seen to the north- 
cast of HD 34078 (A x = -250 - 0", A y = - 300" 
in Figures 1, 4, 5, and 6). This cirrus structure extends 
along filaments roughly from HD 34078 into the region 
containing Filament A**. It may be possible to observe 
this substructure in the optical with high spatial reso- 
lution and a narrow field of view, imaging the nebula 
without strong contamination from the central star (i.e. 
- H ST- ACS). 

These cirrus sheets are the strongest features in the 

4.5 /im IRAC images of this region (Figure 6). The re- 
gion of patchy emission (-330",-75") seen at 8.0 /im is 
observed with a lower relative strength at 4.5 /im. MIPS 
24 /im imaging of this region finds similar relative inten- 
sities as those seen at 4.5 /im, strong cirrus sheets to the 
northeast of the star, with much lower relative intensity 
features coincident with the patchy 8.0 am and "shell" 
structure. This suggests that there is a common dust 
component responsible for the cirrus emission seen at 
4.5, 8.0, and 24 /im, while there is an additional compo- 
nent contributing to the 8.0 /im flux in the "shell" region. 
This additional 8.0 um component is most likely a com- 
bination of an increased contribution from the 7.7 and 

8.6 /im AEFs. One of the results of the IRS scans of the 
bright filaments, discussed below, is that the percentage 
contribution of AEFs to the total flux in a given imaging 
band increases in the brightest knots. We propose that 
the regions showing an enhanced brightness in the 8.0 
/im IRAC band relative to the dust emission at 4.5 and 
24 am are favorable for survival of the AEF carriers. 

Filament A - Filaments A* and A** were covered by 
all four IRAC imaging fields and the MIPS scans at both 
24 and 70 /im. A portion of Filament A was lost at the 
edge of the 70 am map, but covered in all other bands 
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Fig. 8. — Six-color imaging of the Filament A, A*, and A** regions. Data from the four IRAC bands and two MIPS bands are shown 
at their original resolution, with intensity scales intended to display the nebular structure. The black regions on the 3.6 and 70 /im images 
were not observed. The overlays shown on the 8.0 /im image represent the IRS-SL aperture scan direction and the ON/OFF convention 
used in this work. This spectroscopy indicates that the AEFs make an important contribution to the nebular emission at positions of peak 
filament brightness. 



(Figure 8). The 3.6 /im images of this region display a 
peak along the star-facing edge of the clouds containing 
the filaments. The filaments are seen at 4.5 /im, though 
faint compared to the cirrus features seen near the star at 
this wavelength. The 4.5 /im images were the weakest (in 
units of MJy sr -1 ) in our filamentary imaging. We found 
the 4.5, 5.8, and 8.0 /im bands to be tracing the same 
spatial regions, with the 5.8 and 8.0 bands enhanced rela- 
tive to the 4.5 /im. Again, this is attributed to additional 
flux from AEFs in the 5.8 and 8.0 fim bands. 

The 24 and 70 /im MIPS images of Filament A show 
consistency with the 8.0 /im reference image. The 24 /im 
map shows an additional flare structure pointing from 
Filament A** to the northwest (RA = 05 /t 16 m 14.49 s , 
S = +34°24'47.0", J2000). This feature is also seen 
weakly at 4.5 /im. 

Nebular photometry (at 8, 24, and 70 /im) of Filaments 
A and B was performed and the flux densities are pre- 
sented in Tables 2,3, and 4. To ensure a fair comparison, 
the spatial resolution of the 8.0 /im IRAC image and the 
24 /im MIPS scan were degraded. A two-dimensional 
Gaussian kernel was convolved with the data to produce 



8 and 24 /im images with an effective resolut ion of about 
18", comparable to that of the 70 /im map ijRieke et alJ 
I2004|) . The convolved images were integrated over a dis- 
crete number of pixels and multiplied by the solid angle 
subtended by each pixel, 3.4 x 10 -11 sr pixel -1 at 8 /im, 
1.5 x 10 -10 sr pixel -1 at 24 /mi, and 5.8 x 10 -10 si- 
pixel -1 at 70 /im. 

We chose the photometric extraction regions based on 
the sizes of the filaments in the 8.0 /im data. We found 
that the central cores of all four filament positions could 
be contained in a box of roughly 20" x 20". A 20 pixel 
square used at 8 /im, a 10 pixel square was used in the 
24 /im band, and a 5 pixel square extraction was used 
for the 70 /im MIPS data. These extraction regions are 
of the same order as the projected size of the FUSE 
LWRS aperture (30" x 30") used at Filaments A and 
B (Figure 1). The photometric measurements were used 
to constrain the nebular dust temperature in IC 405, as 
discussed in §4.1 and §4.2. 

Filament B - The Filament B region is located 
ps 100 0" north of HD 34078, a part of what iHerbid 
(1958) labels the "E nebulosity". This is a position 
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Fig. 9. — Six-color imaging of the Filament B region. Data from the four IRAC bands and two MIPS bands are shown at their original 
resolution. The images are qualitatively similar, but quantitative analysis shows that the South (HD 34078-facing) edge of filament is more 
extended at 24 /im, suggesting that the carriers of the AEFs are being destroyed in the tenuous intcrfilamcnt medium. The 24 fim map 
also shows a "Plume" to the north and a "Binary" feature to the northeast of the main Filament B emission, not seen as distinctly in other 
bands. Again, the overlays shown on the 8.0 fira image represent the IRS-SL aperture scan direction and the ON/OFF convention. 



TABLE 4 

70 flM MIPS PHOTOMETRY OF NEBULAR FILAMENTS. 



Position 


Program 


RA (2000) 


<5 (2000) 


Flux 






^ h m 


(° ' ") 


(Jy) 


A 


partially off scan map 








A* 


rI6032768 


05 16 19.66 


+34 25 20 


2.562 ± 0.193 


A" 


rI6032768 


05 16 15.64 


+34 23 50 


3.141 ± 0.214 


B 


rI6032768 


05 16 17.74 


+34 32 32 


2.356 ± 0.186 



where FUSE detected fluorescent emission from H2 
as the most prominent spe ctral feature from Ly/3 to 
1190 A l|France et alJ 12004). Filament B was covered 
in all six imaging bands used in this work as well as with 
spectroscopic scans. Unlike the southern filaments de- 
scribed above, all four IRAC bands appear to be tracing 
the same material in the Filament B region. They show 
the same northeast-to-southwest elongated filamentary 
structure with a peanut-shaped central concentration. 
This peanut-shaped central structure, shown in all six 
IR imaging bands in Figure 9, is the center of Filament 
B (RA = 05' l 16 m 16.88 s , S = +34°32'21.0", J2000). 



The MIPS scans at 24 and 70 /im follow the structure 
seen in the IRAC bands, with a few noteworthy excep- 
tions. The MIPS data show that the southern edge of 
the region occurs closer to the star at 24 /im than at 
the IRAC wavelengths. The 24 /im brightness increases 
gradually (B 24 47 - 61 MJy sr" 1 ) over a ~ 30" transition 
region onto Filament B, while the 8.0 /im brightness rises 
more abruptly at the edge of the filament, increasing in 
brightness by more than a factor of two (Bg 18 - 40 MJy 
sr _1 ) over a 15" interface region. Destruction of small 
grains and the carrier of the 7.7 /im AEF in the more 
tenuous environment south of Filament B could produce 
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Fig. 10. — Filament A line strengths measured with the Short- 
Low module on the IRS, as a function of position. The coordinates 
correspond to the 8.0 )im reference image, Figure 1. The map goes 
from East (-245", on the core of the filament) to West (> -220", 
off of the peak of the 8.0 fim filament emission). 
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Fig. 11. — IRS spectra of the ON-filament and OFF-filament 
regions near Filament A. The OFF-Filament was scaled up to the 
ON-filament flux level in the 10 fim region. The ON-filament 
(top panel, blue) spectrum shows the 7.7 and 11.2 AEFs, H2, and 
the 11 — 14 jim "plateau" to be enhanced compared to the OFF- 
filament (top panel, red) pointing. The bottom panel shows the 
difference of the ON and OFF-filament spectra to illustrate the 
differences. 



this effect. We may be seeing the erosion of this cloud 
in progress, the strong UV field of HD 34078 photoe- 
vaporating the carbonaceous grains and weakly shielded 
molecules. This idea could be tested by searching for the 
remnants of these species ([C II] A158 /xni for example; 
van Dishoeck 2004) with future far-IR s pectroscopic in- 
strum entation, such as H erschel-P ACS (Poglitsc h et alJ 
2004). 

The 24 /im MIPS data also reveal a plume of mate- 
rial to the north-northwest of the main knot of Filament 
B (RA = 05H6 m 16.74 s , S = +34°33'09.5", J2000). The 
IRAC images show hints of emission at this location, but 
only at 24 /im is the contrast high enough to identify 
this plume as a separate feature. It is worth noting that 
this plume seems to be "in the shadow" of Filament B, 
analogous to the position of the flare filament seen in 
the shadow of Filament A** at 24 /im. The other no- 
table feature seen in the MIPS imaging is the binary 
structure seen to the northeast of Filament B at 24 /im 
(RA = 05H6 m 21.75 s , S = +34°32'41.9", J2000). The 
two components of this feature are not spatially resolved 
at 70 /im. The "Plume" and "Binary" features are la- 
beled on Figure 8. 

Integrated flux densities for Filament B were measured 
from the IRAC 8.0, MIPS 24, and MIPS 70 /im images. 
The values are listed in Tables 2, 3, and 4. We find the 
24/70 /im ratio to be essentially identical in Filaments 
A*, A**, and B (0.379 ± 0.005), indicative of similar dust 
temperatures in the three regions (§4.2). 

3.2. Spectroscopy: IRS 

Filament A - The IRS scan maps of Filament A have 
been used to measure the variations in emission charac- 
teristics from dense knots to tenuous nebular material. 
The spatial and spectral coverage of the Short-Low (5 - 
15 /xm) module is best suited to studying the change in 



AEF signature while scanning across the region of inter- 
est in full-slitwidth steps. Filament A was mapped from 
east to west, with the eastern portion of the map being 
located on the brightest parts of Filament A, while the 
western portion was off the brightest Filament A emis- 
sion (Figure 10). 

All observed emission lines follow the same generic 
trend: They are brightest near the peak of the 8.0 /xm 
reference image, and their intensity falls off with distance 
from the peak (Figure 10). The AEFs at 6.2, 8.6, and 
11.2 were detected clearly at all Filament A pointings. 
Band substructure in the 6.2 /xm feature is discussed in 
§4.4. Line strengths were measured by integrating the 
spectra over the wavelengths of interest (in the case of 
the 7.7 /im band, the I s * and 2 nd orders were joined 
before measurement). The continuum contribution was 
removed by subtracting a linear fit to the continuum over 
the wavelengths spanned by the line. The 7.7 /im AEF 
was seen to be the most variable spectral feature, rang- 
ing from marginally detectable off the filament to a wide 
(AA ~ 1 /im) band at the peak of the filament bright- 
ness. For all IRS-SL spectral lines, we conservatively as- 
sumed a 20% flux error on each measurement (line and 
continuum) which resulted to total line strength uncer- 
tainties of a few to approximately 30%. 

Figure 11 shows examples of the low-dispersion spec- 
tra acquired at Filament A, showing several aromatic 
and H2 emission lines. The spectra obtained on the peak 
of the 8 /im emission also show a broad "plateau" fea- 
ture from 11-14 /im that is not seen in the spectra 
obtained o ff the filament. T his plateau was first seen 
by IRAS ijCohen et al.1 11985^ . and Spitzer detected a 
similar 10.5 - 14.5 /jm feature in NGC 7023, however 
in that nebula the individual features wash out into 
a br oad AA ~ 2 um fea ture farther from the central 
star ijWerner et all2004b|) . A similar effect with distance 
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Fig. 12. — Filament B line strengths measured with the Short- 
Low module on the IRS, as a function of position. The coordinates 
correspond to the 8.0 )im reference image, Figure 1. The map goes 
from East (-105", off of the filament) to West (> -70", near the 
peak of the 8.0 jim filament emission). 



from the illuminating star was seen by iCesarskv et alJ 
(2000) in the reflection nebula Ced 201. In the case of 
IC 405 however, the entire spectroscopic scan of Fila- 
ment A is roughly equidistant from HD 34078, thus we 
are most likely not observing an evolution with radiation 
field. The imaging indicates that the spectra sample less 
dense regions as the scan moves off the filament. This 
result would be consistent with the 11-14 /im plateau 
being controlled by density 4 . 

Mid-IR spectra of the Orion Bar l)Bregman et all 19891) 
indicate that the molecules responsible for the 11-13 
/im plateau are independent of the carriers of the nar- 
rower AEFs. In IC 405, the AEFs, as well as the lines 
of H 2 and 12.8 /im AEF+[Ne II] (described below), show 
an increase in peak flux with the plateau strength. Fig- 
ure 11 illustrates the changes in these features. The top 
panel shows a comparison between spectra observed at 
the peak of the filament brightness and spectra observed 
off the filament. The "off" spectrum is scaled to the 
10.0 - 10.5 /im region of the "on" spectrum. The differ- 
ence of the two is plotted in the lower panel. The relative 
strengths of the 7.7 /im feature, the H2 lines, and the 11.2 
and 12.8 /an features are all enhanced with the 11 - 14 
/im plateau. 

Several atomic emission lines produced in ionized gas 
were observed in the Short-High observations clustered 
on the brightest region of Filament A, but we only re- 
port on the species that appear in more than 15 of the 
25 pointings. These spectra show strong [Ne II] A12.81 
/im at all positions. In the low dispersion spectra of Fil- 
ament A, the average central wavelength of the 12.8 /jm 
feature is observed at (\ b s ) = 12.817 ± 0.005 /im, con- 

4 We note that while IR, flux density images are representative of 
the underlying density structure, the IRAC flux density does not 
strictly reflect the nebular density jHeitsch e t al.| p006l) . 



sistent with the observed emission coming from ionized 
neon as opposed to the 12.7 /im AEF. In addition to 
[Ne II], [S III] A 18. 71 /mi is detected in most of the high 
dispersion data. We note that it is interesting to find 
these energetic ions (ionization potentials are > 21 eV 
for Ne II and S III) spatially aligned with the emission 
lines of H 2 , described below. One might naively expect 
that the molecules would be destroyed in a region with 
highly excited ions, yet they are observed within the same 
4.7" x 11.3" (IRS Short-High aperture) region. This may 
suggest that the H 2 in IC 405 resi des in dense globular 
reservoirs, as in p l anetary nebulae llHuggins et alJl2002t 
ISneck et all 12003 ILimn et all 1 200ft IHora eiTnTl'lHKi:. 
Narrow band near-IR imaging in the rovibrational emis- 
sion lines of H 2 (A ~ 2 /mi) would be useful in testing 
this suggestion. Finally, we detect emission lines at 12.94 
and 13.14 /im in the majority of the Short-High spectra 
at a signal to noise ratio comparable to the H 2 lines dis- 
cussed below. The 12.94 /im line remains unidentified, 
but we tentatively identify the 13.14 /im feature as the 
O I AA13.139 - 13.157 /im multiplet. Hot pixels in the 
Short-High module do not appear to be the source of the 
unidentified lines, however, the 13.14 /im line does ap- 
pear as a single pixel in some of the pointings. Given 
the possibility of hot pixel contamination and the fact 
that the excitation process for the neutral oxygen line is 
unclear, we do not feel that this is a conclusive detection. 

The pure rotational emission lines of H 2 are seen at the 
on-filament pointings. The (0 - 0) S(5), S(3), and S(2) 
lines are strongest in the core of Filament A, becoming 
weaker in all directions outside of the brightest 20". The 
S(4) line was only weakly observed in one pointing. The 
S(5) line of H 2 is likely blended with [Ar II], as we mea- 
sure a line center that lies between the rest wavelengths 
for the molecular and atomic lines, A fc s = 6.946 ± 0.010, 
where as Ag/gj = 6.91 and X[Arii] — 6-99 /im. The contin- 
uum subtraction for the S(3) and S(2) lines is somewhat 
uncertain due to the variability of the AEFs on whose 
wings they reside. Due to concerns about what effects 
this may have on the determination of the H 2 rotational 
temperature (T(H 2 )), we also used observations of the 
brightest filamentary regions with the Short-High IRS 
module to make clean measurements of the H 2 (0 - 0) 
S(2) A12.28 and S(l) A17.03 /im lines when they are seen 
at high S/N. The determination of T(H 2 ) is discussed in 
§4.2. 

The contributions from the 6.2, 7.7, and 8.6 /im AEFs 
to the 5.8 and 8.0 /im IRAC bands were measured 
in Filament A. Line strengths were determined as de- 
scribed above, and compared to the total spectroscopic 
flux measured in the IR S data over the 5 . 8 and 8.0 /mi 
IRAC imaging windows ijFazio et al.ll2004l I Werner et all 
I2004a|) . We found that the 6.2 /im feature was responsi- 
ble for ~ 33 - 57% of the emission in the 5.8 /im IRAC 
band, with a peak in the AEF contribution coincident 
with the peak in 5.8 /im brightness. The relative con- 
tributions from H 2 S(5)+[Ar II] (A « 6.95 /mi), and the 
7.7 and 8.6 /mi AEFs were measured in the IRAC 8.0 
/im bandpass. The H 2 S(5)+[Ar II] was negligible off the 
filament and rose to a maximum of 3% at the peak 8.0 
/im brightness. The 8.6 /mi AEF contributes 1.0 - 1.5 % 
at all positions, while the 7.7 /im AEF showed a strong 
increase on the filament. At the western end of the IRS 
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Fig. 13. — IRS spectra of Filament B. The blue spectrum in 
the top panel was acquired on the peak of the 8.0 fira emission, 
while the red spectrum was acquired off of the brightest region of 
the filament. The OFF-filament spectrum was scaled to the ON- 
filament flux between 10 - 10.5 fim. The difference of the two 
spectra is plotted in the lower panel. The strengths of the pure 
rotational H2 lines and the 11 — 14 fim plateau are observed to be 
less variable than at Filament A (Figure 11). 



scan (the "OFF" filament position) the 7.7 fim feature 
contributes about 8% of the flux in the 8.0 fim image, 
increasing to 34% at the peak of the 8.0 fim brightness. 
Following this argument and noting the brightness of the 
3.6 fim images relative to 4.5 fim, we consider it likely 
that a substantial fraction of the 3.6 /im brightness is due 
to emission from the 3.3 /im AEF, although we note that 
the scattered starlight component increases at shorter 
wavelengths as well. 

Filament B - The dust, H 2 , and atomic lines observed 
in the Short-Low spectra of Filament B follow two be- 
haviors. The 6.2, 7.7, 8.6, and 11.2 AEFs show a roughly 
linear increase in brightness as the scan moves onto the 
filament, while H 2 S(5)+[Ar II], H 2 S(3) and S(2), and 
AEF+[Nc II] A12.8 fim show a rise in strength then flat- 
ten out across the filament. This is shown in Figure 
12. Examples of the spectra obtained at Filament B are 
shown in Figure 13. The 11 - 14 /im plateau is also ob- 
served at the filament position, although with less vari- 
ability than at Filament A. The Short- High spectra were 
again used to cleanly measure the S(2) and S(l) lines 
near the peak of the Filament B brightness. As in Fila- 
ment A, we detect [Ne II] A12.81 and [S III] A18.71 fim 
at the majority of the pointings. The A12.94 and A13.14 
am features are observed as well. 

Following the procedure outlined for Filament A, we 
measured the contribution of line emission to the 5.8 and 
8.0 fim images. The 6.2 fim AEF is an appreciable frac- 
tion of the 5.8 fim IRAC band brightness, ranging from 
40% at the eastern end of the scan map (OFF filament) to 
46% at the peak of the Filament B. The H 2 S(5)+[Ar II] 
blend makes up 2.8% of the nebular flux off the filament, 
dropping to 1.5% on the center of the filament. This 
is opposite the behavior seen at Filament A, where the 
relative H 2 + [Ar II] contribution increased with filament 
brightness. The 7.7 and 8.6 fim AEFs relative contri- 



bution to the 8.0 fim image increases substantially with 
position. Both lines contribute roughly 2.5 times more 
to the ON filament positions than they do in the more 
diffuse nebular pointings. The 7.7 ^m AEF rises from 13 
to 35% while the 8.6 fim feature increases from 2 to 5%. 

3.3. Regional Comparison 

We can compare the general characteristics of Filament 
B to those of Filament A (Table 5). The H 2 rotational 
emission lines are brighter in Filament A, which could 
indicate a higher H 2 density in Filament A or possibly a 
higher photoelectric grain h eating rate due to the rela- 
tive proximity of HD 34078 ijHollenbach fc TielensllT997l 
lAbel et all 12001 . Similarly, we observe 24 fim emis- 
sion, presumably emission from grains, to be brighter in 
Filament A. Again, we would expect the grain popula- 
tion nearer the illuminating star to display an enhanced 
brightness. However, the AEFs show the opposite behav- 
ior, they are seen to be anti-correlated with the strength 
of the far-UV radiation field. The 6.2, 7.7, 8.6, and 11.2 
fim AEFs all have a position averaged brightness higher 
in Filament B. The 8.0 fim photometric measurements 
are higher in Filament B, consistent with a considerable 
contribution from a brighter 7.7 fim feature. As the AEF 
carriers and small (a < 30 A) grains are thought to be 
heated stochastically by th e absorption of a single UV 
photon (|Draine & Lil l200D. one might expect the AEF 
flux to be higher near HD 34078. We suggest that the 
opposite behavior may be caused by increased photode- 
struction of the AEF carriers in the more intense UV 
radiation field at Filament A. The large column density 
of hydrogen, relative to the aromatics, may place the 
H 2 in a self-shielding regime while the AEF carriers are 
more easily dissociated by the UV radiation field of HD 
34078. An analogous situation is seen in UV illuminated 
PDRs where H 2 is observed while the lower column den- 
sity, unshielded, CO is destroyed l)Luhman fc ,Iaffdll996t 
iLuhman et al]ll997|) . 

IRS mapping observations of the shock region were 
not part of the P20434 program, but predictions about 
the spectral characteristics of the region can be made. 
H 2 (0 - 0) S(3) A9.66 fim, (1 - 0) S(l) A2.12 fim, and 
[O I] A63 fim should be re a dily observable in t his re- 
gion l)van Buren et al.lll990D . iDraine et alJ l) 1983ft make 
predictions about cooling through additional rotational 
lines of H 2 , as well as CO, OH, H 2 0, and [C I]. 

4. DISCUSSION 

4.1. Cometary Bow Shock Associated with HD 34 078 

Far-IR bow shocks are found around roughly 30% of 
galactic O and B runaway stars l)van Buren et al.lll995t 
iNoriega-Crespo et al.ll99"7}l . These stars are ejected from 
the their formation site by interactions, primarily with 
either a supernova binary companion or members of the 
cluste r in which they were formed ijHuthoff & Kaneil 
2002). When one of these high space velocity (> 30 
km s _1 ) stars encounters even moderately dense inter- 
stellar material, a bow shock is formed by ram pres- 
sure confinement of the OB stel lar wind as it moves 
through the ambient medium ijvan Buren fc McCravl 
119881 INoriega-Crespo et al J 11997ft . UV radiation from 
the star heats the shocked gas and dust, producing ex- 
cess far-IR emission in these shock regions, which has 
been observed by IRAS. 
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TABLE 5 

Filament A and B molecular line strength and dust brightnesses. 





Spitzer 
FUSE 


Filament A 
Pos3 


Filament B 
Pos4 


Species 
H 2 




((0 - 0) S(3) Brightness) 3 
((0 - 0) S(2) Brightness) 
UV Brightness 13 
T(H 2 ) (K) c 


2.84 ± 0.62 
1 oo _)_ n no 

.L . _ _ vJ.^O 

2.00 ± 0.60 
418^ 


2.14 ± 0.25 

0.40 ± 0.16 
387if 9 


AEFs 


(6.2 (im Brightness) 
(7.7 /im Brightness) 
(8.6 (im Brightness) 
(11.2 (im Brightness) 


8.73 ± 2.23 
8.52 ± 2.31 
1.69 ± 0.66 
5.05 ± 2.79 


10.13 ± 3.39 
9.39 ± 3.89 
3.30 ± 1.03 
10.58 ± 3.16 


Blends 


(H 2 + [Ar II] 6.95 £tm Brightness) 
(AEF + [Ne II] 12.8 (im Brightness) 


1.23 ± 0.50 
4.04 ± 0.77 


1.29 ± 0.28 
7.35 ± 1.15 


Dust 


F 24 (Jy) 
F 7 o (Jy) 


1.045 ± 0.126 


0.892 ± 0.117 
2.356 ± 0.186 



a H 2 brightnesses are averages over six ON filament pointings at each pos- 
tion. AEF and Blend brightnesses are averages over all IRS-SL pointings 
at each filament. They are in units of (x 10~ 5 ergs cm -2 s -1 sr -1 ). 

b Measured over the ~ 1100 A line complex (1097.7 - 1108.0 A). 

c Measured from the (I(S(3))/I(S(2))) ratio. Filament A shows evidence 
for temperature fluctuations to T(H 2 ) ~ 700 K (see §4.2). 




200 400 
Position (Arcsec) 

Fig. 14. — Dust temperature profiles of the southern portion of 
IC 405, from HD 34078 (at the origin) to Filament A**. The error 
ranges are plotted as dotted lines. The temperature profiles are 
created by taking the ratio of spatial extractions from the 24 and 
70 /im MIPS scans. The dust temperature is seen to rise sharply 
at the location of the bow shock, 10 - 25" north of HD 34078. The 
minimum dust temperature of ~ 60 K is observed at the position 
of Filament A**. This temperature is characteristic of the dense 
filaments in IC 405. 



HD 340 78 was a bow sh o ck ca ndidate in the IRAS 
survey of Ivan Buren et al.| l)1995|) . showing an excess 
brightness at 60 /im, but no shock morphology could 
be detected. A higher spatial resolution ( !') follow up 
to this survey ijNoriega-Crespo et alJl!997j) did not re- 
solve the shock. It is only the imaging capabilities of 
Spitzer (< 6" at A < 24 /im) that allow us to un- 
ambiguously identify the HD 34078 bow shock. Fig- 
ures 4 and 5 show the stellar region at 8 and 24 /im, 



identifying the three shock regions of interest, as de- 
scribed in §3.1. Inspectio n of the images makes clear the 
"cometary" morphology ijvan Buren et al.ll!990(K though 
the ~ 75" diameter arc is more extended than tradi- 
tional ultrac ompact H II regions t hat are interpreted as 
bow shocks ijMac Low et alJll991^ . 

4.1.1. Temperature Structure 

The dust in the bow shock region should be heated by 
the UV radiation field of HD 34078, and we use spatial 
extractions of the 24 and 70 /im MIPS maps to confirm 
this prediction. The F24/F70 ratio cuts have been con- 
verted into dust temperature profiles as a function of 
position following the modified black body prescription 
described in §4.2. The dust temperature profiles for the 
Se, Sc, and Sw regions are shown in Figure 14. These 
cuts extend from roughly 100" south of HD 34078 to 550" 
north of the star, intersecting the Filament A** region. 
We observe the ~ 60 K filament temperature derived 
below for Filament A**, and see the rise in dust tem- 
perature near HD 34078, increasing sharply in the inner 
100" to a peak that is coincident with the location of the 
cometary morphology on the images. The peak dust tem- 

c 

peratures for the shock regions are T d E peak = 82.5±1.9, 

T&o* = 92.3±2.4, and T s d » peak = 91.2±2.5, respectively. 
These peak temperatures are lower limits as the 24 /im 
MIPS data were saturated along a portion of the bright 
bow shock ridge. The dust temperatures derived from 
the photometric measurements (given in Table 1) of the 
shock region average over the peak 24 /im brightness and 
the surrounding emission. 

The dust temperatures we derive are factors of 2 - 3 
high er than seen in the origina l IRAS bow-shock sur- 
vey ijvan Buren fc McCravlll 988). The high temperature 
of the HD 34078 shock can be related to the high velocity 
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Fig. 15.— Short-High spectra of the H 2 (0 - 0) S(2) and S(l) 
pure rotational emission lines at 12.28 and 17.03 fim. An offset of 
0.4 Jy was added to the Filament B spectra for this figure. Ratios 
of the measured line strengths (I (S (2)) / 1 (S (1))) are used to derive 
the H 2 rotational temperature, T(H 2 ). The 7(S(2))//(S(1)) tem- 
perature in Filament B agrees with that derived from Short-Low 
measurement of the I (S (3)) / 1 {S (2)) temperature (T S (H 2 ) sd 400 
K). However, the Short-High observations of Filament A show ev- 
idence for non-uniform excitation corresponding to a temperature 
of roughly 700 K. 



of the star and the high spatial resolution of these obser- 
vations. The ~100 km s _1 velocity of HP 34078 is out- 
side th e range of velocities considered bv lvan Buren et alJ 
(1995). HD 34078 has the largest one-dimensional proper 
motion (within the errors) of these shock candidates 
(prriy in Table 2 of van Buren et al. 1995). The high 
spatial resolution provided by Spitzer-MlPS allows us 
to spatially resolve the hottest shock regions. The larger 
IRAS beam (5' for van Buren et al 1995; and 1' for 
Noriega-Crespo et al. 1997) sampled the entire shock re- 
gion, diluting the warmer inner regions with cooler dust 
on the exterior of the shock. Averaged over 5' around 
HD 34078 (-90 - 210" of Figure 14), the temperature 
derived for Sc decreases to 76 K. 

4.1.2. Bow Shock Density: H2 Excitation and UV Dust 

Extinction 

The morphological and thermal evidence places the HD 
34078/IC 405 interface in the cometary bow shock fam- 
ily, and we can use the analytic models of van Buren & 
McCray (1988) and van Buren et al. (1990) to explore 
other topics related to the IC 405 system. An estimate 
of the nebular density enables us to calculate the column 
density and UV optical depth associated with the bow 
shock and explore the implications for existing studies. 

The distance from the star to the shock region can be 
calculated from Equation (2) of (van Buren et al. 1990; 
see also van Buren & McCray 1988), assuming an exist- 
ing knowledge of the nebular density. Alternatively, we 
can measure the star-shock distance directly from the 
IRAC and MIPS observations, and solve for the density. 
Let I be the measured distance from the star to the shock 
structure in the direction of motion, «15" on the 8.0 /im 
IRAC image. At the assumed distance of 450 pc, this 
corresponds to a separation of 1.77 x 10 15 cm. Rewrit- 



ing equation (2) of Ivan Buren et all l)1990j) . the density 
{riH = 2n# 2 + uhi + njjn) in the region of IC 405 
through which HD 34078 is currently moving is given by 



n H 



(3.03 x 10 38 ) m»,_ 6 v WtS ^ v~l r z cm -3 (1) 



where to*^ is the stellar wind mass loss rate in units of 
x 10~ 6 M Q yr _1 , v Wt g is the terminal velocity of the stel- 
lar wind in units of x 10 s cm s _1 , u*^ is the stellar veloc- 
ity in units of x 10 6 cm s _1 , and p,n is the dimensionless 
mass per H-nucleus. We take hh t o be 1 and assume 
a ste l lar velocity of ~ 1 00 km s _1 ( Bl aauw &: Morgan! 
Il95l iBoisse et al.l l2005h . The stellar mass loss rate 
(log(m») = -9.5) a nd terminal wind ve locity (800 km 
s _1 ) are taken from lMartins et alJ l)2005l) . Following this 
procedure, we calculate a nebular density of tih — 244 

cm~ 3 . 

I i ■ v ' : Equation (7) of lvan Buren et aTl l)199flh the 
column density of the swept-up material is given by 



N 



(1.95 x 



inl9\ • 1/2 1/2 -1/2 -1 1/2 -2 

10 m ' „ v' A Urr v, fi riri cm 



6 u w.8 r"H u *,6 H 

(2) 

with the parameters defined as before. This calcula- 
tion yields a column density of material piled-up by the 
intrusion of HD 34078, N = 4.8 x 10 17 cm" 2 . This 
va lue is very s i milar to the H2 column density found 
bv lBoisse et alJ l)2005|) for the hot absorption component 
along the HD 34078 line of sight. They find N(H 2 , J > 
5) = 3.3 x 10 17 cm~ 2 , comparable to their model pre- 
diction of 2.4 x 10 17 cm -2 . Our derived column density 
of the bow shock region supports the hypothesis that the 
high-lying states of H2 observed on the HD 34078 sight- 
line are produced in a hot layer near the star, swept-up 
by the passage of HD 34078, and excited by stellar UV 
photons. 

The sightline to HD 34078 is seen to be hi ghly reddened 
comp ared with other regions of IC 405 l)France et alJ 
2004). This is inferred from the blue rise of the ratio 
of nebular surface brightness to stellar flux (S/F») at 
far-UV wavelengths (900 - 1400 A). The Spitzer data 
allow us to determine if this enhanced reddening can be 
attributed to the ma terial associated with the sho ck re- 
gion. Equation (6) of|van Buren fc McCravl (JT988) gives 
the UV optical depth {tjjv) of the swept-up dust. As- 
suming that the majority of the dust survives the inter- 
action with HD 34078 and that the swept-up material 
has a velocity of a few km s _1 , Tjjy is of the order 0.01. 
This relatively small value of Tjjy is insufficient to ex- 
plain the observed S/F* in IC 405. In other words, the 
dust responsible for the differential extinction towards 
IC 405 is not associated with the bow shock. It is more 
likely that this dust is located in the translucent cloud 
along the HD 34078 sightline. This cloud is observed to 
produce the major ity of the cold (J = 0, 1; T(H 2 ) = 77 
K) H 2 absorption ( Bois se et al]l2005^ . and we conclude 
that this cloud is also responsible for the majority of the 
extinction towards HD 34078. 

4.2. H2 and Dust Temperatures in Filaments A and B 

Observationally, little attention has been paid to IC 
405 in comparison with other bright PDRs (NGC 2023, 
NGC 7023, IC 63, and the Orion Bar). Here we establish 
the dust and H 2 rotational temperatures of the outer 
filaments. 
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H2 - Rotational temperatures can be measured 
from H2 emission line ratios in the IRS band- 
pass (jRosenthareFal] HqH \Me£setj^fXMk ■ This 
temperature directly probes the bulk of the molecular 
mass in PDRs and has been a widely used diagnostic 
since a large number of mid-IR observations were made 
available by ISO (see references in the review by van 
Dishoeck 2004 and §1.1). We have used the IRS Short- 
Low measurements of the position averaged H 2 S(3) and 
S(2) emission line strengths and Short-High observations 
of the central cores of Filaments A and B to determine 
T(H 2 ) for these regions. Only the brightest six Short- 
Low pointings corresponding to the ON position in Fig- 
ures 10 and 12 were considered. The use of the Short- 
High module allows us to separate the (0 - 0) S(2) A12.28 
/jm and S(l) A17.03 /im lines from nearby spectral fea- 
tures. The S(2) and S(l) line strengths were measured 
from the two cleanest (defined as a combination of line 
S/N and a low, featureless background) Short-High ob- 
servations at each filament. The S(2) and S(l) spectra 
from Filaments A and B are shown in Figure 15. 

T(H2) was determined by combining the measured line 
ratios, I(S(2))/I(S(l)) and I (S '(3))/ 1 (S '(2)) , with en- 
ergy levels and transition probabilities from the litera- 
ture. This temperature determination assumes that the 
pure rotational lines are optically thin and thermally 
populated, with an equilibrium ortho-to-para H2 ratio of 
3. Internal and foreground attenuation is also assumed 
to be negligible, which is supported by the finding that 
the reddening towar ds the nebula is es sentially zero at 
far-UV wavelengths llFrance e t alJl20 04D. 

The energy levels from iMandv fc Martini l)1993() give 
AE2i/k = 667 K, and the Einstein- A radiative tran- 
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Fig. 16. — FUSE observations of Filaments A (above) and B 
(below). H2 fluorescence is observed at each position imposed on 
the dust scattered continuum of HD 34078, providing direct evi- 
dence for a population of photo-excited nebular gas. A synthetic 
spectrum of H2 emission, scaled to the flux level of the data, is 
overplotted. Lyman band H2 (labeled — and 1 — 0) and photo- 
spheric absorption lines are also observed. 



sition probabilities are A g( 2 ) _ 2.76 x 10 9 s 1 and 
A s(1) = 4.76 x 10" 10 s" 1 llWolniewicz et al.lll99l . The 
ratio of level degeneracies, gs(2)/9sa), is 9/21. Using 
these parameters, we obtain rotational temperatures of 
T^(H 2 ) = 692 ± 55 K and Tf x (H 2 ) = 381 ± 28 K for 
Filaments A and B, respectively. These temperatures 
are higher than one would expect for a cold interstellar 
cloud, but are similar to recent PDR observations that 
find typical H2 rotational temperatures in the range of 
a few hundred degrees K ijThi et alJH999t IHabart etall 
l20T)llAllers et al.ll2005l) . 

Noting the added uncertainties in the continuum lev- 
els discussed above, we also calculated the H2 rota- 
tional temperature using position averaged strengths of 
the S(3) and S(2) lines observed in the low-dispersion 
spectra. Including the values for ^5(3), AE^/k, 
from the references given above with gs(3), we found 
T3VH2) = 418±73 18 K and Tf 2 (H 2 ) = S87±H K. The 
agreement at Filament B is remarkable (381 vs. 387), 
but the Filament A temperatures differ by almost 300 
K. 

We interpret variations seen at Filament A as differing 
relative populations in adjacent rotational states, indica- 
tive of the influence of non-thermal excitation by far- 
UV photons. The distribution of excited robvibrational 
states produced in the fluorescent cascade differs from a 
thermal distributio n and this is reflected in the o bserved 
emission line ratios ( Black & van Dis hoeckHl987|) . If the 
filament structure is clumpy, the narrower field-of-view of 
the Short-High aperture (as opposed to the position av- 
eraged Short-Low data used to determine T32(H2)) could 
be dominated by a knot with a larger non-thermal emis- 
sion fraction. Two tests of the Short-Low data have been 
performed to explore this possibility. We have used the 
individual H 2 S(3) and S(2) line strengths to look at T 32 
as a function of position. This shows that the derived 
H 2 rotational temperature displays excursions to greater 
than 700 K, albeit at the lower signal to noise positions 
not considered in the average line strengths used to de- 
termine T32 above. Secondly, we looked at temperature 
variations along the Short-Low slit at a pointing that 
spatially overlapped the Short-High spectra used. We 
measured a higher temperature (T 3 4 2 (H2) > 520 K) along 
the star-facing edge of Filament A, coincident with the 
Short-High positions, suggesting that the line ratios used 
to determine the H2 rotational temperature may be al- 
tered by UV photons in this region. 

The majority of rovibrational emission lines that would 
indicate fluorescent pumping fall at wavelengths shorter 
than the Spitzer-IRS bandpass, thus IRS spectroscopy 
alone may not be sufficient to constrain the excitation 
mechanism of nebular H 2 . In IC 405 however, direct ev- 
idence for non-ther mal (fluorescent) exc itation of H2 has 
been presented bv iFrance et all l)2004D , FUSE obser- 
vations clearly revealed the far-UV spectrum of fluores- 
cent molecular hydrogen at Filaments A and B. These 
spectra are shown in Figure 16. Filament A is shown 
in the top panel, where the emission from H2 is seen 
imposed upon dust-scattered light from HD 34078. We 
note that the fluorescent lines at Pos3 display a broaden- 
ing in excess of the instrumental width, as has been seen 
in th e H 2 emission spectra of other PDRs ijFrance et alJ 
I2005|) . HD 34078 is clearly identified as the excita- 



16 



France et al. 



tion source by the highly excited photospheric features 
of P V and Si IV seen in the scattered light. As has 
been shown for fluorescent H2 em ission near the Trapez- 
ium ({France fc McCan dlissll2 005l) . a star with such pho- 
tospheric ions l|Pellerin et alll2002T ) is the most likely to 
have the requisite far-UV flux (A < 1110 A) to produce 
observable UV fluorescence. The FUSE spectrum of Fil- 
ament B is shown in the bottom panel, over plotted with 
a synthetic H2 spectrum l)France et al.ll2004j) . It is inter- 
esting to note that the synthetic spectrum was created 
with a rotational temperature of 400 K, only attempting 
to find qualitative agreement with the data. The fluo- 
rescent H2 was integrated over the strongest far-UV line 
complex (1097.7 - 1108.0 A) after subtracting the con- 
tribution from dust scattered light. Integrated surface 
brightnesses are quoted in Table 5. 

Fluorescent excitation by HD 34078 can also be ob- 
served in the highly excited H2 absorp tion lines on the 
HD 34078 sightline l|Boisse et al.ll2005li . As non-thermal 
emission clearly plays a role in IC 405, we propose that 
the average rotational temperatures of the filaments is 
~ 400 K, while UV-pumping has produced non-uniform 
excitation. The result of this excitation has been to al- 
ter the observed line ratios at the Short-High pointings 
used in Filament A, leading to the observed temperature 
deviation. 

Dust - Thermal dust temperature (T^) can be deter- 
mined from mid an d far-IR photometry made w idely 
availa ble by IRAS l|Sodroski et alJ fl987t iSchnee et al.l 
2005). Measurements of T^ fit two or more IR data 
points, assuming emission from the dust cloud is de- 
scribed by a modified black-body curve. The black-body 
function is typically modified by a power-law emissivity 
function that is representative of the grain composition. 

We have used the 24 and 70 MIPS observations 
of IC 405 to derive the dust temperature in the three 
filaments where the scans overlap (A*, A**, and B), as 
well as across the bow shock region (§4.1.1). The con- 
tribution of discrete spectral features to the photometric 
flux in these bands has been considered. The weak spec- 
tral features in IRS Long-Low (20 - 30 /im) observations 
contribute less than ~ 8% of the 24 /im band flux at all 
pointings. The contribution of line emission to the 70 
/xm brightness is more difficult to constrain without the 
aid of spectroscopy. We note that [O I] A63 /im lies at 
the edge o f the MIPS 70 /im fi lter. [O I] may be strong 
in IC 405 lHvan Dishoeckll2004l). but without direct mea- 
surement of the nebular SED at these wavelengths, we 
do not attempt to correct for emission lines. 

The observed F24/F70 photometric ratio, computed 
from Tables 3 and 4, was fit using a modified black 
body spectrum, with a dust emissivity proportional to 
v 2 . The choice for the dust emissivity law was deter- 
mined by the theoretical dust absorption opacity wave- 
length dependence 2.92xl0 5 (A//iin)~ 2 , valid in the 
20 ix m< A <700 /im range for the average diffuse 
ISM <|Li Draind 120011) . The spectral slope in the 24 
and 70 /im filter bandpasses does not make a significant 
contribution to the resulting brightness ratios, the cor- 
rection being significantly smaller than our measurement 
errors for T d > 40 K. 

The F24/F70 ratios were nearly identical in the three 
filaments, 0.379 (§3.1), and the inferred dust tempera- 
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Fig. 17. — IRS measurements of the ratio of 6.2 and 11.2 /im 
line strengths in Filaments A and B. Their ratio is an indication of 
the ionization state of the aromatics in these regions. The relative 
ionization fraction increases off the core of Filament A, but is seen 
to be constant at Filament B. Note : The Filament B data is 
displayed ON (-50") to OFF (-105") for comparison with Filament 
A in this plot. 



tures reflect this. The mean temperature for the neb- 
ular dust in IC 405 is 63 K. Individually, the ob- 
served filament dust temperatures are T^ = 62.8 ± 1.5, 
Tf* = 62.6 ± 1.4, and Tf = 62.7 ± 1.6 K. These 
grain temperatures are roughly a factor of 3.5 higher than 
those found for molecular clouds in Perseus, Ophiucus, 
and Serpens {(T d ) = 17 K; Schnee et al. 2005), sug- 
gesting that the grains have processed the UV photons 
from HD 34078. Our derived temperature for the IC 405 
filaments is similar to the lower limit of the dust tem - 
perature (55 K) in the Orion Bar PDR ijLis et al.lll998|) . 
What may be surprising about the IC 405 dust measure- 
ments is the similarity of the temperature at Filaments 
A* and A** (roughly equidistant from HD 34078) with 
that of Filament B. Single photon heating effects should 
be important for grains with < 1000 carbon atoms in th e 
intense radiation field of HD 34078 iJDraine inil200l . 
This may be the cause of the constant temperature dis- 
tribution in the nebular filaments. 

4.3. Ionization and Molecular Environment 

There is general consensus that the AEFs are produced 
by a superposition of vibrational modes of l arge carbona- 
ceous molecules l|Verstraete et al.l [2001; Pcctcrs ~et all 
2002). Laboratory and theoretical studies have shown 
that the typical interstellar aromatic emission spectrum 
is composed of a combination of neutral and ionized 
carriers with a distribution of sizes (see the review by 
van Dishoeck, 2004, and references therein). These 
studies have shown that cations (positively charged 
ions) have enhanced 6-10 /im features arising from 
stronger C - C stretching modes. Neutral molecules 
have stronger C - H bands, producing the major- 
ity of the aromatic emission at 3.3 and 10 - 14 
/im (jHonv et al.l200lUBregman k, Temil2005D . The ratio 
of 6.2/1 1.2 /im emission bands can be indica tors of ion- 
ization <)Uchida et all2000UBakes et al.l200H) and molec- 
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Fig. 18. — The ionization measure of the AEF carriers (6.2/11.2 
/jm) as a function of strength of the molecular hydrogen emission 
at 9.66 (im. The correlation observed at Filament A suggests that 
ionization of the aromatics increases in lower density environments 
at given radiation field I Bakes ct al. 2001). No trend is observed 
at Filament B. 



ular structure <)Chan et alJl200H iHonv et alJl200H) . 

IRS scans of the two filaments in IC 405 allow us 
to study the change in ionization state with environ- 
ment. Assuming that Filaments A and B are located 
at roughly the same heliocentric distance as HD 34078, 
we can constrain the relative radiation fields at the two 
positions. Variations in the spectral properties along 
Filaments A and B can be correlated with change in 
filament structure/density and changing radiation field. 
Previous studies have primarily focused on the variations 
of AEFs with radiation f ield strength llVerstraete et alJ 
Il996t IPeeters et al.ll2003 iWerner et aljl2004b|) Our re- 
sults for molecular ionization fraction as a function of 
position are shown in Figure 17. The scan moves off the 
filament to the right, the position coordinates have been 
calculated to correspond to the images presented in this 
work. The 6.2/11.2 ratio is displayed as squares for Fil- 
ament A and triangles for Filament B. We observed an 
increase in ionization fraction as the scan moves from the 
central knot to the diffuse nebula at Filament A, going 
from roughly 1.4 to 3.3. This can be understood as an 
increase in the contribution from ionized ca rriers in less 
dense regions at a constant radiation field ijBakes et alJ 
2001). 

The ionization fraction is seen to be constant with po- 
sition in Filament B. The baseline 6.2/11.2 /mi ratio is 
lower in Filament B than in Filament A (~ 0.9 vs 1.4), 
suggesting that the overall ionization fraction is smaller 
in Filament B. Filament B is roughly twice the angular 
distance from HD 34078 and consequently receives about 
a quarter of the UV flux present at Filament A, con- 
sistent with the FUSE observations presented in §4.2. 
While the ionization fraction declines with diminishing 
UV flux, there might not be a one-to-one relationship. 
This conclusion is derived in connection to the assertion 
made in §3.3 that Filament A shows indirect evidence 



for enhanced AEF carrier destruction relative to Fila- 
ment B. Using the Filament B region as a baseline, our 
observations are consistent with the photodestruction of 
the aromatic molecules responsible for the AEFs being a 
multi-step process where dissociation is preceded by ion- 
ization. Additional observations of PDRs are available 
in the Spitzer and ISO archives, and we will explore 
similarities with IC 405 in future work. 

The AEF ionization measure can also be correlated 
with the environment in which H2 emission is produced. 
A comparison of the H2 emission strength and AEF ion- 
ization is shown in Figure 18. The S(3) rotational line 
was used as representative of H2 emission because of its 
brightness and separation from other lines. We observe 
the AEF ionization fraction to increase as the molecular 
hydrogen line strength decreases in Filament A. Assum- 
ing that H2 is representative of density, these results also 
indicate that the AEF carrier population is more ion- 
ized in more tenuous r egions at a given UV radiation 
field IjBakes et alj|20 01). Combining this result with the 
spatial dependence of the molecular ionization described 
above confirm that in strongly illuminated PDRs, the 
6.2/11.2 fim ratio is a good diagnostic of the local ion- 
ization cond i tions. This result is consistent with the 
iBakes et al.l l|200lD prediction for regions with a ratio 
of ionizing radiation field to electron density (G /n e ) of 
~ 10 4 cm 3 . No correlation is observed at Filament B. 

4.4. Substructure in the 6.2 fim Feature - Evidence for 

PANHs? 

In addition to providing information about the ioniza- 
tion state, the 6.2 /mi fe ature can be a diagn ostic of the 
AEF carrier structure. IPeeters et alJ l)2002j) have per- 
formed a survey of the 6.2 /mi feature in a variety of as- 
tronomical objects, and have determined empirical sub- 
classes for this band. Lines with peak wavelengths near 
6.22 /mi form Class A, while lines peaking between 6.24 
and 6.3 /mi make up Classes B and C. Most H II regions, 
reflection nebulae, and non- isolated Herbig Ae/Be stars 
belong to Class A, while isolated Herbig Ae/Be stars and 
most planetary nebulae belong to Class B. They find red- 
peaked (A ~ 6.3 /mi) lines to be composed of "pure" 
PAHs, while blue-peaked (A ~ 6.2 /mi) lines have N- 
a toms substituted at in ner lattice sites. 

IHunVins et all {2005) carry the nitrogenation hypoth- 
esis further, determining the necessary site for N- 
substitution in order to reproduce the observed 6.2 /tm 
wavelength and ruling out increasing molecular size alone 
as the explanation for the observed line profiles. They 
find that nitrogen substitution is required to reproduced 
the observed profiles, labeling this PAH subclass poly- 
cyclic aromatic nitrogen heterocycles (PANHs). They 
conclude that PAHs or PANHs with N-substitution on 
the edges of the molecule produce the emission at A > 6.3 
/mi, while PANHs with the N-atom substituted in the 
interior of the mole cule are responsible for the Class A 
spectra observed bv IPeeters et all ((2002) . 

With this observational and theoretical evidence in 
mind, we report a tentative detection of band substruc- 
ture to the red of the main 6.2 /tm peak at Filament A. 
Filament A spectra are shown in Figure 19 from 5.8 - 
7.2 /tm. The spectra are dominated by a main Class A 
feature at X max — 6.210 ± 0.007 /tm . At the peak fil- 
ament brightness position, we observe an enhanced red 
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FIG. 19. — IRS Short-Low observations of the 6.2 fim AEF 
in Filament A. The spectra have been normalized to the peak of 
the on-filament 6.2 /an band and offset by factors of 0.5 for this 
display. The data show a variation in the 6.2 ^m profile from the 
core of the filament (ON, top) to the ed ges (OFF, bottom). The 
spectra maintain a Class A 6.2 /im peak IPeeters et al.H2002fl . but 
develop red substructure away from the core of the filament. Three 
Gaussians are fitted to the spectra showing the main peak (blue), 
the second peak seen at 6.33 fim (green), and the broad 6.5 (im 
"hump" (red). The 6.5 Gaussian on the second OFF spectrum 
(offset = +0.5) was fit over a slightly expanded wavelength range 
to show the feature more clearly. 



wing to the main feature. This wing develops a distinct 
second peak at 6.33 /jm as the overall 6.2 band decreases 
in strength away from the central core of the filament. 
An additional red "hump" appears further to the red 
(A ~ 6.50 fim) becoming stronger with increasing dis- 
tance from the peak filament brightness position. Three 
Gaussians are shown overplotted on the data in Figure 
18 representing each component. The H2+[Ar II] blend 
is seen near 6.95 fim. 

The red features are qualitatively similar to the 
coronene cation (C 2 /i H^~r,) spectra shown in Figure 3 of 
IHudgins et al.l l|2005|) . The substructure is most pro- 
nounced in the same region where the ionization fraction 
is seen to increase. One interpretation is that process- 
ing by UV photons is altering the PAH population in 



the interface between the dense filament and the more 
tenuous nebular regions. We only see red substructure 
in one Filament B spectrum. Filament B is generally 
observed to have an unmodified Class A 6.2 fim band, 
A mra = 6.219 ± 0.004 fim. We note that a 6.35 fim fea- 
ture was reported as a red shoul der of the 6.2 um ban d 
in ISO observations of M17-SW (jVerstraete et alJl996|) . 
interpreted as large PAHs un der the influence of a strong 
UV radiation field. However. IPeeters et al.ll|2002D report 
no detection of a 6.35 fim feature. 

5. SUMMARY 

We have presented new Spitzer observation of IC 405, 
an emission/reflection nebula in Auriga illuminated by 
the runaway 09.5 star HD 34078. We have used IRS 
spectroscopy to assess the morphological differences seen 
in four band (3.6, 4.5, 5.8, and 8.0 fim) imaging obtained 
with IRAC. The combination of the IRAC data with 
24 and 70 fim MIPS maps of the nebula reveal wave- 
length dependent morphology in several subregions. We 
attribute these differences to variations in the relative 
strengths of AEFs and emission from grains. In addition 
to the AEFs, the spectroscopy revealed emission from 
the pure rotational lines of H2 and atomic fine structure 
lines. The emission bands/lines were observed to vary 
in both absolute and relative strength as a function of 
position at two bright nebular filaments. 

Dust temperatures in three nebular filaments were de- 
rived from 24 and 70 fim photometry. The dust temper- 
ature was a constant 63 K, independent of the distance 
from the illuminating star. A similar temperature anal- 
ysis performed on spatial profiles of the 24 and 70 fim 
data obtained near HD 34078 show elevated dust tem- 
peratures (> 90 K) in a bow shock ahead of the star, in 
the direction of stellar motion. The bow shock scenario is 
supported by the cometary morphological structure seen 
in the 8.0 and 24 fim images of the stellar region. We 
use this data to support the two-component model for 
the HD 34078 sightline. The excited H 2 is thought to be 
associated with the shock region, while the majority of 
the cold H2 absorption and dust extinction occurs in an 
intervening translucent cloud. 

The measured AEF line strength were found to be anti- 
correlated with the intensity of the far-UV radiation field 
in the nebula, possibly due to increased destruction of 
the aromatics. The spectroscopic data also were used to 
constrain the relative ionization state of the AEF car- 
rier molecules. The ionization state was seen to increase 
with decreasing density in one filament. This ionization 
fraction was seen to be anticorrelated with the H2 emis- 
sion line strength, consistent with enhanced AEF carrier 
ionization with decreasing density at a given UV radia- 
tion field. Evidence was presented for substructure in the 
6.2 fim AEF, suggestive of changing carrier size and/or 
composition as a function of nebular environment. Line 
ratios were used to derive the H2 rotational temperature 
in two filaments, indicating an average rotational tem- 
perature of roughly 400 K, with evidence for additional 
non-uniform excitation in the filament nearer to the illu- 
minating star. These measurements were complemented 
by observations of the far-UV H2 fluorescence spectrum, 
providing direct evidence for excitation by UV photons. 

As new imaging and spectroscopic observations of 
PDRs become available in the Spitzer archive, oppor- 
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tunities for comparative studies with existing UV data 
are made possible. Combining the UV and mid-IR ob- 
servations with near-IR observations of the rovibrational 
emission lines of H 2 , an accounting of the energy redistri- 
bution following electronic excitation can be carried out 
in order to test models of H 2 fluorescence. A panchro- 
matic H2 spectrum could be produced in a variety of 
PDRs, from pure reflection nebulae (NGC 2023) to star- 
forming regions (Orion Bar) to H II regions with harder 
radiation fields (planetary nebulae such as M 27) . These 
data would allow models that span a wide range in ex- 
citation temperature, density, and radiation field to be 
compared with observation. H2 observations can be com- 
pared to the AEF and dust properties in these regions 
in order to explore the relationship between molecular- 
phase PDR diagnostics and local environmental charac- 
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